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Abstract 



Several lines of evidence point to a relationship between gamma-ray bursts (GRBs) 
and the high mass stars that explode as supernovae. Arguments that GRB sources 
accelerate cosmic rays (CRs) are summarized. High-energy neutrino detection from 
individual GRBs would mean that they are hadronically dominated, that is, that the 
Q\ \ amount of energy deposited in the form of nonthermal hadrons in GRB blast waves 

greatly exceeds the radiated energy inferred from their photon emission spectra. 
Such a detection would make GRBs the favored candidate sources of ultra-high 
energy and super-GZK CRs. Cascade radiation induced by high-energy hadrons 



could be detected from GRBs with 7-ray telescopes. 
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1 Introduction 



Considerable evidence connects long-duration GRBs to star-forming regions 
(see (U I2) for review) and, consequently, to the high-mass stars that radiate 
strongly in the UV band and evolve to supernova core collapse. For example, 
the host galaxies of GRBs have blue colors, consistent with galaxy types that 
are undergoing active star formation. GRB counterparts are found within the 
optical radii and central regions of the host galaxies. Lack of optical counter- 
parts in approximately one-half of GRBs with well-localized X-ray afterglows 
(so-called "dark bursts") could be due to extreme reddening from large quan- 
tities of gas and dust in star-forming host galaxies. Supernova-like emissions 
have been detected in the late-time optical decay curves of a few GRBs in- 
cluding, most recently, the low redshift (z = 0.17) GRB 030329 @). 

Furthermore, X-ray features have been detected in about 6 GRBs, which im- 
plies the existence of dense, highly enriched material in the vicinity of the 
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sources of GRBs (see (]4j; |5|) for review). This includes the variable absorption 
feature attributed to photoionization of Fe in GRB 990705, Ka and recombi- 
nation edges from highly ionized Fe and S in the afterglow spectra of GRB 
991216, multiple high-ionization emission features detected in GRB 011211 
and, at lower significance, claimed X-ray Fe Ka features in GRB 970508, 
GRB 970828, and GRB 000214. 

The association of GRB 980425 and the Type Ic SN 1998bw 0), if true, di- 
rectly connects GRBs and SNe. The well-measured light curve of SN 1998bw 
has been used in a number of cases to model reddened excesses in the opti- 
cal afterglow spectra found tens of days after the GRB, again pointing to a 
SN/GRB connection. 

Because supernovae are widely thought to accelerate galactic cosmic rays, the 
supernovae associated with GRBs should play a special role in the CR/GRB 
connection. Here we summarize work connecting CRs, SNe, and GRBs, be- 
ginning with source models of GRBs and arguments for CR acceleration by 
supernovae. High-energy neutrino detection with km-scale telescopes such as 
IceCube would conclusively establish that hadrons are accelerated by GRBs, 
and is likely if GRB sources are hadronically dominated by a factor of ~ 10 
or more. This would be required in a GRB model that accelerates nonther- 
mal power-law distributions of particles to produce ultra-high energy cosmic 
rays (UHECRs) with energies £ 10 18,5 eV, and super-GZK CRs with energies 
£ 10 20 eV. Hadronic cascade radiation could make GRBs luminous sources 
of GeV-TeV emission that would be detectable with GLAST, VERITAS, and 
HESS. 



2 Gamma-Ray Bursts and Supernovae 



Leading scenarios for the sources of long-duration GRBs are the collapsar and 
supranova models. In the former case (|7j), GRBs are formed in the seconds to 
minutes following the collapse of the core of a massive star to a black hole. 
During the collapse process, a nuclear-density, several Solar-mass accretion 
disk forms and accretes on the newly born black hole at the rate of ~ 0.1- 
1 M Q s _1 to drive a baryon-dilute, relativistic outflow through the surrounding 
stellar envelope. The duration of the accretion episode corresponds to the 
period of activity of the relativistic winds that are argued to produce the 
prompt variable GRB emission. A wide variety of collapsar models are possible 
(0), but their central feature is the one-step collapse of the core of a massive 
star to a black hole. A major difficulty of this model is to drive a baryon-dilute, 
relativistic outflow through the stellar envelope (0). 

In the supranova model floT). GRBs are the second in a two-step collapse 
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process of an evolved massive stellar core to a black hole through the interme- 
diate formation of a neutron star with mass exceeding several Solar masses. 
The neutron star is initially stabilized against collapse by rotation, but the loss 
of angular momentum support through magnetic dipole and gravitational ra- 
diation leads to collapse of the neutron star to a black hole after some weeks to 
years. A two-step collapse process means that the neutron star is surrounded 
by a SN shell of enriched material at distances of ~ 10 15 -10 17 cm from the 
central source. The earlier SN could yield ~ 0.1-1 M of Fe in the surround- 
ing vicinity to provide a surrounding shell of material that forms prompt and 
afterglow X-ray features in GRB spectra. 

A pulsar wind and pulsar wind bubble consisting of a quasi-uniform, low den- 
sity, highly magnetized pair-enriched medium within the SNR shell is formed 
by a highly magnetized neutron star during the period of activity preceding 



its collapse to a black hole (11). The interaction of the pulsar wind with the 
shell material will fragment and accelerate the SNR shell, and the pulsar wind 
emission will be a source of ambient radiation that can be Comptonized to 
gamma-ray energies (jl2l) . The heating of the SN shell by the plerionic emission 
could give rise to a characteristic cooling shell signature, which might explain 
some instances of the delayed reddened excesses in optical afterglow spectra 
0). 

In both the collapsar and supranova models, the stellar progenitors are £ 1OM 
stars that evolve to core collapse and produce a supernova with an expanding 
nonrelativistic shell, and a GRB with its highly relativistic ejecta. The shocks 
formed in these nonrelativstic and relativistic outflows can accelerate cosmic 
rays. 



3 Cosmic Rays and Supernovae 



The strongest argument that hadronic cosmic rays (CRs) are powered by SNRs 
is probably the claim that only SNe inject sufficient power into the Galaxy to 
provide the measured energy density of CRs |l3 : 0). A time- and volume- 
averaged kinetic luminosity ^ 10 41 ergs s _1 in the disk of the Milky Way is 
required to power the galactic CRs. The galactic SN luminosity is ~ (1 SN/30 
yrs) x 10 51 ergs/SN w 10 42 ergs s _1 which, given a 10% efficiency for converting 
the directed kinetic energy of SNe into CRs that seems feasible through the 
shock-Fermi mechanism, is adequate to power the cosmic radiation. 

Besides available power, Fermi acceleration at SNR shocks offers another 
strong argument that cosmic rays originate from supernovae. In Fermi models, 
the directed kinetic energy of a nonrelativistic or relativistic outflow is trans- 
ferred to a small fraction of suprathermal particles, either through repeated 
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scattering and diffusion upstream and downstream of the shock front (first- 
order Fermi acceleration), or through energy diffusion accompanying stochas- 
tic pitch-angle scatterings with the magnetic turbulence spectrum (second- 
order Fermi acceleration). 

As compared with some other acceleration mechanisms, for example, charge- 
depletion fields in pulsar magnetospheres or through magnetic reconnection in 
stellar flares, the shock-Fermi mechanism naturally produces power-law par- 
ticle spectra, at least in the test particle approximation. First-order Fermi 
acceleration at a nonrelativistic strong shock gives the canonical test-particle 
shock index s = —2, and first-order Fermi at a relativistic shock gives the 



canonical test particle shock index s = —2.2 (16). Modifying these standard 
injection indices through energy-dependent diffusion for impulsive, steady, or 
stochastic injection models offers a wide range of flexibility to explain obser- 
vations of cosmic rays. 



4 Cosmic Rays and Gamma-Ray Bursts 



Arguments for a connection between CRs and GRB sources are based on the 
coincidence between available power of GRBs within the GZK radius and 
the power required to accelerate super-GZK particles. Moreover, the available 
power from GRBs within our Galaxy is sufficient to accelerate CRs between 
the knee at « 3 x 10 15 eV and the ankle at ~ 10 18 ' 5 eV (if these CRs are trapped 
in the halo of our Galaxy), and could even power a significant fraction of the 
GeV-TeV/nuc cosmic rays f)l7l ). 

The Larmor radius of a particle with energy 10 20 -E 2 o eV is ~ 100E 2 o/ (ZB^ G ) 
kpc. Unless super-GZK CRs are heavy nuclei, they probably originate from 
outside our Galaxy. Based on assumptions that the faintest BATSE GRBs 
are at redshifts z « 1, it was shown (Q El) that the energy density of super- 



GZK CRs is comparable to the energy density that would be produced by GRB 
sources within the GZK radius, assuming that the measured 7-ray energy from 
GRBs is roughly equal to the total energy deposited in the form of UHECRs. 
This coincidence is verified by a detailed estimate of GRB power in the context 
of the external shock model for GRBs . A prediction of this hypothesis 

is that star-forming galaxies which host GRB activity will be surrounded by 
neutron-decay halos. 

The local density of L* galaxies like the Milky Way galaxy can be derived from 
the Schechter luminosity function, and is ~ 1/(200-500 Mpc 3 ). The BATSE 
observations imply ~ 2 GRBs/day over the full sky. Due to beaming, this rate 
is increased by a factor of 500 (|2lh. This implies an event rate of about 1 GRB 
in the Milky Way every 10 3 -10 4 years, and an energy injection rate of ~ 10 40±1 
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ergs s -1 (22). This estimate takes into account the assumption that the GRB 
rate density follows the star-formation rate history of the universe, and a factor 
of ~ 3 to account for the contribution from clean and dirty fireballs, such as 
the X-ray rich GRBs, for every classical long-duration GRB. Thus the GRB 
rate is about 10% as frequent as Type Ib/c SNe, and about 1% as frequent as 
Type II SNe. 

Because each GRB has a total energy release of a few xlO 51 ergs, a factor 
of a few greater than the energy release in "normal" SNe, but occurs ~ 100 
times less frequently, the available power from the sources of GRBs is a few 
per cent of the power from Type la and Type II SNe. The efficiency for accel- 
erating CRs in the relativistic outflows of GRBs could be much greater than 
in nonrelativistic SNe, so that the progenitor sources of GRBs could make a 
significant contribution to CRs at all energies. 

From the rate estimates, we see that about 1 in 10 to 1 in 100 SNRs would 
exhibit this enhanced emission from strong CR acceleration due to an associ- 
ated GRB. The better imaging and sensitivity of the GLAST telescope and 
the next generation of imaging ground-based air Cherenkov telescopes, namely 
VERITAS, HESS and MAGIC, will test this hypothesis. 



5 High-Energy Neutrinos from GRBs 



Detection of high-energy neutrinos from km-scale telescopes such as IceCube 
will establish hadronic acceleration in sources such as blazars and GRBs. We 
(J23I ) have recently completed the first quantitative set of calculations of pho- 
tomeson neutrino production for the collapsar and supranova models of GRBs 
that also treats associated 77 pair-production attenuation of the emergent 
7-ray emission. This model takes into account nonthermal proton injection 
followed by photomeson energy loss, based on our photo-hadronic model for 
blazar jets (J2J). The basic assumption in our calculations is that equal ki- 
netic energy is injected in nonthermal protons as is measured in hard X-ray 
and soft 7-ray photons. Given the fluence of a GRB, the energy injected in 
protons therefore depends only on the Doppler factor 5. 

Neutrinos will be formed through photomeson production with internal syn- 
chrotron radiation and, in the supranova model, with external synchrotron 
photons emitted by the pulsar wind electrons. Even with the very strong 
radiation field that could occur in the supranova model, model-independent 
arguments (jiih . which we reproduce here, show that it is not possible under 
the assumption of equal energy in hadrons and observed radiation to detect 
neutrinos from any GRBs except those exceedingly rare events at the level 
£ 3 x 10~ 4 ergs cm -2 , which occur about 2-5 times per year (j25f). 
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Fig. 1. Energy fluence of photomeson muon neutrinos for a model GRB (23). The 
thin curves show collapsar model results where only the internal synchrotron radi- 
ation field provides a source of target photons. The thick curves show the 5 = 100 
and 300 results for the supranova model calculation, which includes the effects of 
an external pulsar wind radiation field. 

The detection efficiency in water or ice of ultrarelativistic upward-going muon 
neutrinos (i/ M ) with energies e u = 10 14 er eV is P V(l = 10~ 4 e^., where \ — 1 



for €t < 1, and x = 0.5 for > 1 (26). For a neutrino fluence spectrum 
parameterized by v<§> v = 10~ 4 

with a km-scale v detector such as IceCube with area A v 
therefore 



4 ey" erg cm 2 , the number of v ^ detected 



10 10 A 10 cm 2 is 



NJ> e T ) 



d* ^ - 0.6 ^ 



g av (e T ) , for e T < 1 
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(1) 



for e T £ 1 , 



where g av { e T) = 1 + [(1/2°^) — 1](1 — e^"). For a v<& v spectrum with a v ~ 0, 
the number of to be expected are N v ~ 1.20_ 4 yiio(l + ^lney 1 ) for er < 1, 
and iVy ~ 1.2^_4ylio/ v / eT for ey > 1. Thus, if the nonthermal proton energy 
injected in the proper frame is comparable to the radiated energy required 
to form GRBs with hard X-ray/soft 7-ray fluences £ 10 -4 ergs cm -2 , then 
extremely bright GRBs are required to leave any prospect for detecting 
with km-scale neutrino detectors. This estimate does not take into account the 
loss of efficiency to produce from energetic protons through photomeson 
processes, so the number of detected from individual GRBs would be even 
less. 

Fig. 1 shows the total fluences expected from a model GRB with a photon 
fluence of 3 x 10~ 5 ergs cm -2 that is radiated in 50 one-second pulses. The vari- 
ability time scale determines the size scale of the emitting region given 8, and 
thus the intensity of internal synchrotron photons for photomeson production. 
The thin curves show collapsar model results at 8 = 100, 200, and 300. The 
expected numbers of that a km-scale neutrino detector would detect are 
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N v = 3.2 x 10~ 3 , 1.5 x 10~ 4 , and 1.9 x 1CT 5 , respectively. The heavy solid and 
dashed curves in Fig. 1 give the supranova model predictions of N u = 0.009 
for both 5 = 100 and 5 = 300. The equipartition magnetic fields used in this 
calculation are 1.9 kG and 0.25 kG, respectively. The external radiation field 
in the supranova model makes the neutrino detection rate insensitive to the 
value of 5 (as well as to the variability time scale, as verified by calculations). 
As the above estimate shows, there is no prospect to detect from GRBs at 
this levels or from a GRB with a factor of ten more fluence. 



Although this calculation is not optimistic for neutrino detection, note that 
the coincidence between the power released by GRBs in 7 rays within the 
GZK radius and the power required for the super-GZK CRs assumes that 
equal energy is injected into the highest energy cosmic rays as is detected as 
GRB photon radiation. This is possible in a second-order Fermi acceleration 
scheme where a very hard, very high energy particle population is produced 
(271). If particles are instead injected with a —2 or —2.2 spectrum through a 
first-order Fermi mechanism, then GRBs must be hadronically dominated if 
they are to power the UHECRs. Thus, ~ 10-100 times as much power must 
be injected in hadrons as is inferred from the hard X-ray /soft 7-ray emission 
from GRBs. Under these conditions, IceCube could expect to detect high- 
energy neutrinos from several GRBs each year. Such a discovery would have a 
profound impact on our understanding of cosmic rays. In addition, luminous 
MeV-GeV-TeV cascade radiation induced by hadronic secondaries would be 
produced. Observations of distinct 7-ray components in GRB spectra would 
also support the GRB/CR connection. 



6 Conclusions 



Many lines of evidence point to a close linkage between supernovae, gamma-ray 
bursts, and cosmic rays. Detection of high-energy neutrinos from the sources of 
GRBs would confirm hadron acceleration to very high energies. Under the as- 
sumption that equal energy is injected into protons as is measured in radiation, 
we have found (23) that there is no realistic prospect for neutrino detection 
except from the most fluent GRBs at a level ^3 x 10~ 4 ergs cm -2 . More 
optimistic estimates from the viewpoint of detecting GRB neutrinos could be 
found in proton-dominated GRB models. Without this hypothesis only the 
brightest GRBs can be expected to be detected with both high-energy 7-ray 
and neutrino detectors. The assumption of hadronically dominated GRBs is 
required to provide the required energy in UHECRs if the protons are acceler- 
ated in the form of a power law with number index steeper than —2. Detection 
of high-energy neutrinos from GRBs would test scenarios for GRBs, the hy- 
pothesis that GRBs sources are powerful accelerators of ultrarelativistic CRs, 
and our understanding of the origin of GRB radiation. 
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